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Abstract 

— -J ■ Ion Beam Sputtering (IBS) is known to produce surface nanopatterns over macroscopic areas 

Q , on a wide range of materials. However, in spite of the technological potential of this route to 

o ■ 

nanostructuring, the physical process by which these surfaces self-organize remains poorly under- 
K^ stood. We have performed detailed experiments of IBS on Si substrates that validate dynamical 

\^ • and morphological predictions from a hydrodynamic description of the phenomenon. Our results 



elucidate flow of a nanoscopically thin and highly viscous surface layer, driven by the stress created 
by the ion-beam, as a description of the system. This type of slow relaxation is akin to flow of 
macroscopic solids like glaciers or lead pipes, that is driven by defect dynamics. 
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I. INTRODUCTION 

In the paper that historically coined the word "Nanotechnology" [ij, erosion of solid 
targets through ion-beam sputtering (IBS) was already put forward as the most promising 
technique to structure the surface of a wide range of materials. Indeed, for ion energies 
ranging from 100 eV to 100 keV, the procedure has shown a remarkable capability |2| to 
produce ordered nano-scale sized patterns (mostly ripples and dots) over large areas (up to 
;ens of cm^) for a wide range of targets, including semiconductors, metals, and insulators 
3|, y]. However, despite its large potential for technological applications J5[, the promise 
of IBS as a fully controlled method to tailor patterns with custom designed properties has 
turned elusive. This is not only due to the challenges that description of the dynamics of 
surfaces in times of order 1 s pose to multi-scale modeling 6| when microscopic events take 
place in times of order 1 ps. Rather, it originates in the lack of a basic understanding of the 
physical nature of the mechanism that controls the pattern formation. 

Classically, since the seminal work by Bradley and Harper (BH) [7| the interplay between 
sputtering and surface diffusion had been identified as the key mechanism leading to pat- 



tern formation in IBS. Thus, a characteristic length scale would be selected 



from the 



competition between the morphologically unstable dependence of the sputtering yield with 



Iq 



local surface curvature, and thermal surface diffusion that smooths out surface features |9|. 
However, thus far only partial qualitative agreement had been reached between this classical 
description and experiments |3|. Only after the recent realization of the non-trivial role of 



impurities in the emergence of the pattern 



10 



ll| has the need arisen to dispose of the BH 



mechanism for the simplest case of monoelemental semiconductor targets like silicon, that 
are amorphized by the ion beam 12|, as in recent experiments on Ar+ irradiation |l3|, Il4l ]. 
For instance, one of the most direct implications of the BH picture, that ripple formation 
should take place for any incidence angle 9 between the ion beam and the normal to the 
uneroded target, has been experimentally disproved, there being a critical angle 9c such 
that ripples appear only for 9 > 9c. Incidentally, for ion energies between 10 and 40 keV, 
higher than in [l3|, |lj], this fact was already noted by Carter and Vishnyakov (CV) 15| . 
who explained it phenomenologically via a smoothing effect of momentum transfer from the 
ions to the target atoms. 

To date, two main non-BH-type mechanisms have been proposed to explain formation 



or absence of patterns in IBS of monoelemental semiconductors [l6|. Namely, mass re- 
distribution g, ll7H19| and ion induced solid flow 20|. The former employs results from 
Molecular Dynamics (MD) simulations in order to rephrase the CV effect as the influence 
of material (rather than momentum) displacement due to the beam, on the surface mor- 
phology. The latter proceeds, rather, through a continuum description of the surface flow 
that is driven by the surface confined stress due to the accumulated damage produced by 
irradiation. Remarkably, both descriptions agree with experiments with respect to e.g. the 
value of 9c and the dependence of the pattern wavelength on 9 16|. Note that, as far as 



physical mechanisms are concerned, MD simulations [6|, |17H19| are fundamentally limited 
in the IBS context, since they cannot probe the required macroscopic time scales at which 
non-trivial surface evolution takes place. This limitation has required e.g. the ad-hoc use 
of surface confined viscous flow 21| in p, [l8|, in order to obtain the correct time scales for 
the description of the dynamics. What remains to be done is, rather, to prove if a physical 
description that incorporates viscous flow from the outset is indeed able to predict properties 
of the pattern that can be experimentally tested. 

In this paper, we perform IBS experiments of Si in order to validate specific predictions 
from the hydrodynamic description put forward in 20|, that stem from the assumption that 
the main driving field in the process is the ion-induced residual stress that builds up in 
the flowing layer. Consistency between experiments and theoretical predictions thus allows 
us to elucidate the physical mechanism through which pattern formation occurs in these 
non-equilibrium nanoscopic systems. 



II. PREDICTIONS OF THE THEORY 



The main idea of the solid flow description of IBS 20[ is that, as a consequence of the 
impact of the ions and the subsequent release of energy within the target, defects are created 
inside the material. These events occur in a few picoseconds after the impact. Relaxation 
of some defects leads to sputtering of target atoms, but also to the generation of a residual 
stress that is confined to a thin amorphous layer that builds up beneath the surface 22] and 
reaches a stationary thickness. This ion-induced (compressive) stress is characterized by a 
slow time relaxation which involves a highly viscous flow of the amorphous layer, that will 
be assumed to be incompressible, in a nanometric analogue of the motion of glaciers or lead 
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FIG. 1. Schematic view of an IBS experiment. An incident energetic ion (red) impinges onto the 
target at an angle 9, inducing a coUision cascade that amorphizes a thin region through creation 
of vacancies and interstitials. At long time scales, the amorphous solid flows as a highly viscous 
fluid. Here 7 is the local slope of the surface profile h{x, t) (see the main text for notation and 
further details about the equations). 



)ipes 23| . The effect of these slow events can be cast into a hydrodynamical description 
24j | of the form (see Fig. [T]for a schematic description of the system): 



V ■ V = 0, V ■ T = Vp, 

n • T ■ n = -CTK + Ti^^*\ 



n ■ T ■ t = V(7 ■ t + T, 



(ext) 



V(x,t) 



0, 



(1) 

(2) 
(3) 
(4) 



where V and p are fluid velocity and pressure, T is the stress tensor inside the amorphous 
layer, n (t) is the unitary vector locally normal (tangent) to the interface, a is the interface 
surface tension and k its curvature. Equations ([2]) an d (p) are evaluated at the surface and 
Eq. (jlj) at the amorphous-crystalline (ac) interface 20|]. In our range of interest for ion 
energies (say, 300 to 1000 eV), the ripple wavelength (tens of nm) is much larger than the 
thickness of the amorphous layer (a few nm) 3|, |j] . Introducing this fact into the analysis of 
Eqs. ([ID-dl]) amounts to taking the so-called shallow-water lirait 2J]. Using further a linear 



approximation in perturbations around a flat target profile, we obtain a real part for the 
linear dispersion relation (namely, the rate at which periodic perturbations with wave vector 
q grow or decay |8[) that is given by 

where (j){9) = dg {^{0) sin(6')), d is the average thickness of the amorphous layer, and fi is 
its (ion-induced) viscosity. The parameter /e can be understood as the gradient of residual 
stress induced by the ions across the amorphous layer, whose angular dependence is described 
through the function \l/(^). Since this angular function needs to be prescribed, we take 
^(0) = cos 9 as the simplest geometrically motivated choice that shows good agreement [20|] 
with previous experiments, as well as with the ones presented here. 

Given \E' and j'e, we can extract testable predictions from the theory through the value 
of the ripple wavelength. A, occuring to linear approximation of Eqs. ([l])-(jl]). 



A = 2vr,/-^. (6) 



The first prediction has to do with the angular dependence of the patterns. For ^1^(6*) = cos 6 
[hence, (f){6) = cos(26')], we find a value for the critical angle 6c = 45°. Moreover, this value 



does not depend on the ion energy, in good agreement with X-ray experiments 18 



and in 



contrast with an energy- dependent critical angle from MD simulations in 6|, from 6c — 40° 



at 250 eV to 6c — 35° at 100 eV 25|]. Note, experiments at low energies {e.g., 100 eV) are 
difficult to perform reliably with a small dispersion. 

Thus far, the analytical description is purely kinematic, in the sense that we have not 
made explicit which is the physical mechanism behind the stress. Experimentally, local 
measurements of stress in a nanoscopic amorphous layer is problematic if not unfeasible. 
Thus, a natural procedure is to detail the nature of stress generation in the theoretical 
description, relating it with relevant physical parameters of the experiment, such as average 



ion energy, E, or flux, J, and obtain testable predictions from such an assumption. In [26|] a 
simple (experimentally validated) model was proposed based on knock-on implantation and 
stress relaxation by defect migration to the surface. This leads to a dependence of the stress 
that builds up throughout the amorphous layer, T oc fsd, on ion energy of the form 



Here, Y and z/ are the material Young's and Poisson's moduli, and R, p are other material- 
dependent constants (related to properties like the binding energy). From Eq. ([7]) it follows 
that, for low energies, T ~ JE^^'^, while for large energies (which depend on R, but typically 
E > 100 eV), T ~ -E"^/^. We also need a prescription for the dependence of the amorphous 
layer thickness on energy. It is customary to assume a power law dependence of the form d ~ 
E"^"^, where m is assumed to be in the range 1/3 — 1/2 and can be determined approximately 



27j . Thus, for a fixed angle, /e oc T/d ~ E '^1^ ^™, and using 



from TRIM simulations 

we find 

A ~ E", (8) 

with a G [0.92, 1.08], namely, the dependence of the ripple wavelength with energy is roughly 
linear. 

Another prediction can be obtained from Eq. ([5]), that allows to find the characteristic 
scale for the exponential growth of the pattern amplitude occurring at short times during 
with linear approximation holds. Thus, 

r{E, 9) ^ 2vr/c.,.,. ~ -^^^^^ ~ -^^^, (9) 



with h G [1.33, 1.67], q* = 27r/A, and we have used that /i ~ 1/J (2l|. Notice that e.g. the 
surface roughness, W, grows in the linear unstable regime as H^ oc e*/^. This exponential 
behavior will eventually be interrupted by nonlinear mechanisms (coming, for instance, from 



stress 20| or from purely erosive effects 28|) at sufficiently long times. Hence, the time 



duration of the validity of the linear approximation is also controlled by r. 

III. EXPERIMENTS 

Si (100) targets (380 m thick, p-type, 1 — 10 Ohm cm) were sputtered with Ar"*" ions 
under different incidence angles within the 300 — 1000 eV energy range. The angle 6 was 
adjusted by axial rotation of the sample in front of the ion gun with an overall resolution 
of ±2°. The ions were extracted from a commercial 3 cm beam-diameter Kaufman-type ion 
gun (VEECO) located 25 cm away from the target. In order to avoid metal contamination 
on the surface during irradiation, the sample holder was covered with a sacrificial Si wafer 
and the sample was attached to it with a double sided conductive vacuum tape. Prior to 
the process, the current density at the sample position in the plane parallel to the source 

6 



grids {9 = 0) was set to 50 /iA/cm^ with a Faraday cup located on a movable shutter that 
prevents, additionally, unwanted irradiation of the sample. The surface morphology of the 
irradiated surfaces was imaged ex-situ with a Nanoscope Ilia equipment (Bruker) operating 
in intermittent contact mode and using silicon cantilevers (Bruker) with a nominal radius 
of curvature of 10 nm. 

In Fig. [2] we show four AFM micrographs of Si (100) surfaces bombarded a.t E = 700 eV 
and different incidence angles. The surface does not contain any visible pattern for 6 = 45° 
(or for smaller angles, not shown). On the other hand, as seen in Fig. [3^, roughness values 
increase significantly above this critical angle |l4l |. 




(c) 



(d) 



FIG. 2. Experimental morphologies of ion irradiation of Ar"^ on Si. All the figures are obtained at 
constant dose, at energy E = 700 eV, and at angles of incidence, 6 = (a) 45°, (b) 55°, (c) 65°, (d) 
75°. Note the absence of pattern at 45°. Lateral size is 2 fiui for all the panels, with ions arriving 
from right to left in all cases. 



In order to compare experimental data with the model predictions derived above, we 
note that the linear regime is expected to have a duration that changes with experimental 
conditions, see Eq. (Q. For instance, it is expected to last shorter for increasing angles 
of incidence 6 > 6c (in a form reminiscent of critical slowing down for continuous phase 
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FIG. 3. (a) Roughness dependence on incidence angle for experiments performed at constant dose 
(empty diamonds) and at the intrinsic time scale predicted by Eq. ([9]) and estimated through 
(jlOP (solid black circles) for E = 700 eV. (b) Time evolution of the surface roughness (circles) for 
E = 700 eV and 9 = 55°. The blue solid line is an exponential fit (linear regime), while the green 
dashed line is a power law fit with exponent 0.32. (c) Time evolution of the ripple wavelength 
(squares) for the same conditions as in (b). The time regime in which the wavelength grows with 
time (coarsening) starts at t ~ 25 minutes as indicated by the dotted line; this is another signature 
of the onset of non-linear effects. At shorter times the wavelength is almost constant. The same 
crossover time separates exponential from power law behavior in (b). 

transitions |8|, [l6|). Actually, Eq. (Q allows us to control the experimental times, in order 
to guarantee that the system is truly evolving within linear regime, given a fixed reference 
experiment in which such state can be unambiguously assessed. Thus, given a pair of angle- 
energy reference values, (6'ref, -E'rcf), we can extrapolate the value of r for any other pair 
(9, E) through 

^(Q r^\-^(0 p ^ ■/cxp(^rcf,^rcf)^ref' 0^(^rcf) ..,„s 

T{t>, E) - T{Oreh ^^refj ; .. ^x ^ 7_^o^ ,,,,, , (lUj 



J, 



Gxp ^^7 



E)E' 



'l+2mjJ2 



(0) 



where Jexp(6',-E') = Jcxp{E)cos9 is the flux used in a particular experiment at energy E 
and angle 9, with J^^piE) being the flux at normal incidence. In our case, we choose 
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Oy-cf = 65° and E'ref = 700 eV. Thus, from Eq. (fTOj) we can extract the experimental times 
for different angles or energies which we define as intrinsic time scale, namely, the time at 
which experiments at different angles and/or energies are within linear regime. 

In Fig. [3] we show how the surface roughness depends on the conditions under which 
the experiments have been performed. Thus, we emphasize that doing the experiments at 
constant time (or fiuence) may produce surfaces which are described by different regimes 
(linear or non-linear), depending on the value of 6. The morphologies in Fig. [2] correspond 
to the times obtained at constant dose in Fig. [3] where typical nonlinear motifs, like facets. 



are recognized, special 



y in Figs. |2t-d. Note how the onset of these effects correlates for 



different observables [^ (power law growth of W, coarsening of A) in Figs. |3)d,c. 
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FIG. 4. (a) Ripple wavelength dependence on the angle of incidence for E = 700 eV. The empty 
diamonds correspond to perpendicular ripples like those in Fig. [2]3, obtained at constant dose for 
all angles. The empty triangle is for parallel ripples with the same energy. Full symbols correspond 
to experiments performed in the proper linear regime as estimated by the intrinsic timescale in 
Eq. Q. The solid blue line is a fit to Eq. ^, the green dashed lines corresponding to the same 
fit but assuming a ±2° uncertainty in 6. The full black circles correspond to perpendicular ripples 
and the full red squares to parallel ripples, (b) Ripple wavelength dependence on the ion energy 
for constant dose experiments (empty diamonds) and intrinsic timescale experiments (solid black 
circles). The solid lines are fits to a power law with exponent 0.47 and a linear dependence (as 
noted in the figure). 

In addition, in Fig. H^ we show the dependence of the ripple wavelength on the incidence 
angle. The empty symbols stand for experiments performed at constant dose and the full 
ones are for experiments at the intrinsic timescale defined in Eq. (j9]) that guarantees linear 
behavior. The blue solid line corresponds to a fit to Eq. (Q using \&(^) = cos 6*, and the 
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green dashed lines are a similar fit taking into account a ±2° uncertainty in the experimental 



measurement of 9 29|. These fits confirm the validity of Eq. ([6]). Moreover, in Fig. HJd we 
show the dependence of the wavelength on the ion energy. Again, the prediction of the 
theory that A scales almost linearly with ion energy fits nicely, provided the experiments are 
performed at the intrinsic timescale, and not at (the customary) constant dose experiments. 
In general, agreement is not reached for the latter since those data points correspond to the 
nonlinear regime that sets in for times longer than the scale ([9]), at which predictions made 
from linear approximation break down. This fact proves the self-consistency of our present 
analysis and stresses the predictive power of the solid flow theory. 

IV. DISCUSSION AND CONCLUSIONS 

We have presented theoretical predictions based on ion induced residual stress and viscous 
flow as the main physical mechanism driving surface nanopattern formation by IBS, together 
with the experimental validation of such predictions. These results allow us to extract 
important consequences. Note that the predictions in terms of ion energy originate from 



the physical model in 26|, in which direct knock-on implantation is assumed to produce 
stress that can be relaxed by defect migration to the free surface of the amorphous layer. 
When taken as an input of our hydrodynamic framework, this allows us to provide the 
scaling of different observables with energy. Thus, the ripple wavelength is predicted to scale 
roughly linearly with E^ in marked contrast with the classical E'^'"^ scaling from BH theory 
7|. Other theoretical descriptions of the problem, both based on effective interface 







equations (see [16| for an overview) or on atomistic simulations combined with multi-scale 
approximations, are based on mere superposition of mechanisms that actually operate at 
different time scales, rather than on a single framework that leads to a dynamic description 
as is the case here. As an example, in the original BH theory sputtering takes place in a 
few ps, while surface diffusion occurs in ns. One could actually expect contributions from 



sputtering to surface-confined transport [2l|, |28|, but adding them together ad-hoc, as is 



natural to a first approximation, will necessarily overlook them. 

Finally, an important part of our analysis relates to the dependence of the linear or 
nonlinear behavior with experimental parameters and observation time. Thus, an essential 
requirement for experimental reproducibility and meaningful comparison with theory is to 
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make sure that measurements made for different conditions (e.g. energy, angle, etc.) corre- 
spond to the proper dynamical regime. We expect that clarification of the physical basis of 
IBS allows to enhance experimental control over the technique, that finally brings it up to 



fl 



the well-founded high expectations expressed almost 40 years ago [l| 
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